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DEFINITION OF SYMBOLS 




Symbol Definition 

t temperature (°F) 

T temperature (*R) 


pressure 



(denoted as a function of temperature) 


P' 


T 


dP T 

~dT 


; also used simply as p' 


~ dT 


j ; also used simply as p" 


/ Btu 1 

specific Gibbs' function; ie, ^ = h- Tsl 




specific enthalpy 


Btu 

lbm 


specific entropy 


vp ^ two-phase specific Gibbs' function; (denoted as a function of 

temperature) 


d* r 


dT 


; also used simply as 


d 2 ^» 

^"t dT^ ’ a * so usec * s ^ m piy as 

v specific volume (ft 3 /lbm) 

u specific internal energy (Btu/lbm) 
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DEFINITION OF SYMBOLS (Cont'd) 


Symbol 



cr 

v 

a 

a 

M 

V 


Definition 


values of above functions on saturation liquid curve; 
all are denoted as temperature functions 


change of above functions after complete isothermal 
phase change; all are denoted as temperature functions 


single-phase constant volume specific heat capacity (Btu/lbm 0 F) 

single-phase polytropic specific heat capacity 

two-phase constant volume specific heat capacity 

single-phase sonic velocity (ft/sec) 

two-phase pseudo-sonic velocity (ft/sec) 

pseudo- Mach number (V/a) 

fluid flow velocity (ft/sec) 
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DEFINITION OF SYMBOLS (Cont’d) 


Symbol 

q 

A 

G 

W 

x 

L 

t 3 

r 

f(T) 


f'(T) 


J 


a 


b 

A 

T 


Definition 

specific heat transfer ( Btu/lbm) 
cross-sectional area of duct (ft 2 ) 

mass flow rate per unit cross-sectional area (lbm/sec-ft 2 ) 
total mass in system 
quality of two-phase mixture 
flow system length (ft) 

absolute triple point temperature (°R) (491. 69 for H 2 0) 

ratio of absolute temperature in ®R to triple point temperature 

an arbitrary function of temperature formed from the partial 
integration of the internal energy function derivative 

df(T) 

dT 

ft/lb f 

conversion factor 778. 26 — — — 

Btu 

, . _ „„ lbm ft 

conversion factor 32. 17 — — r 

lb^ sec 

an empirical coefficient (489. 286 Btu/lbm for H 2 0) 
an empirical coefficient (489. 281 Btu/lbm for H 2 0) 
polytropic exponent ( 1. 252 for II 2 0) 
time (seconds) 


viii 
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DEFINITION OF SYMBOLS (Concluded) 

Symbol Definition 

Subscript 

f denotes liquid phase 

g denotes gas or vapor phase 

t denotes throat conditions 

1 denotes a reference state; arbitrarily taken as the duct inlet 

state in numerical examples 

Superscript 

* denotes choke conditions for Fanno or Rayleigh flow 
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PSEUDO-SONIC VELOCITY AND PSEUDO-MACH NUMBER 
CONCEPTS IN TWO-PHASE FLOWS 

SUMMARY 


In ideal gas theory , the introduction of the Mach number into the 
analysis of steady, one-dimensional duct flow results in associating choking 
conditions with a Mach number of unity for several well known flows , such as 
Rayleigh and Fanno pipe flows and isentropic converging-diverging duct flow. 
Replacing the ideal gas by a liquid-vapor, two phase mixture in the above 
named flows complicates the analysis even when equilibrium between the phases 
is assumed. However, it is possible to introduce a pseudo-sonic velocity 
concept with a correspondingly defined pseudo-Mach number concept such that 
these three two-phase flows exhibit choking conditions at a pseudo-Mach number 
of unity. 


It is further shown that a point-by-point, non-tabular analysis of 
the above two-phase flows is possible when the equations for the Gibbs' function 
and the vapor pressure curve are known. In the case of the liquid-vapor region 
of H 2 0 an empirical equation for the Gibbs' function is derived, while the vapor 
pressure curve is developed theoretically by the imposition of a transformation 
property normally encountered in ideal gas theory. 


INTRODUCTION 


The use of the pseudo-sonic velocity will be delayed until a few 
thermodynamic relations can first be reviewed and specialized for two-phase 
usage prior to their incorporation into the definition of the new variable. 


For a closed system, the differential of the Gibbs' function defined 


by 


* = h - Ts (1) 

is given by 

d'l' = ~ vdp - sdT . (2) 

tJ 

In this form the pressure and the temperature are the independent variables 
except when the pressure depends only upon the temperature, as is the case in 
equilibrium, two-phase thermodynamics. For this case if the pressure 
dependence upon the temperature is expressed by 

P = P T 

then 

dp - p'rpdT (4) 

where p' T is the slope of the vapor pressure curve at the temperature T. If 

dp from this equation is substituted into equation (2) the temperature derivative 
of the Gibbs' function becomes 

*' T = J vp' T - s . (5) 

That the right side is indeed a temperature function becomes apparent when 
the two-phase equations 


v 



+ 



s 



+ XV 


f Sn 


and the Clapeyron equation 



(6) 

(7) 


(8) 


2 



are substituted into the right side of equation ( 5) to give 




T 




(9) 


When equation (5) is used to define the entropy in a two-phase region 
as a function of T and v, the entropy differential becomes 


ds 




+ - p* dv 
J v T 


Hence for a reversible process 


dq 



+ - Tp'dv . 
J 


(10) 


(ID 


From this equation it is seen that the coefficient of dT represents a theoretical 
two-phase specific heat capacity at constant volume. It is designated by the 
symbol cr^ to distinguish it from the usual single-phase c ^ symbol. In general, 

the chart of two-phase properties representing the T-s plane shows that a 
curve of constant volume possesses discontinuous slopes at the saturation 
envelope. This means that the saturation value of cr^ is not the same as the 

saturation value of c^. Thus, for any point inside the two-phase region 


<7 

V 




(12) 


As a final digression before defining the pseudo-sonic velocity it 
is recalled that the velocity of propagation of a small disturbance in a non- 
condensable single phase is given by 



where ol is the acoustic or sonic velocity. In terms of the specific volume equation 
( 13) becomes 


a 


2 



(14) 


3 


V 


By direct analogy, the right side of equation (14) will be defined as 
the pseudo-sonic velocity parameter for two-phase analysis and the symbol 
”a" will be used for its description ; thus 



To evaluate the partial derivative in this relation, equations (4) and (12) are 
used to transform equation ( 10) into 


v j 

ds = rrv dp + r p'dv . (16) 

lp' J 

From this the desired partial derivative is readily obtained and substituted into 
equation ( 15) to give for the square of the pseudo-acoustic velocity 


a 2 = 


j Tv 2 p' 2 

|c 

J O’ 


(17) 


These relationships are discussed further in Appendix C. 

For positive a the right side is likewise positive; hence the pseudo- 

sonic velocity parameter is real when cr^ is positive, which is the usual case. 

This may be readily checked by observing the slopes of constant volume curves 
in the T-s plane. As examples, these slopes are positive for the liquid-vapor 
regions of water, nitrogen, para and normal hydrogen, etc. 

From equations ( 17) and ( 12) it is apparent that "a" depends ex- 
plicitly upon T and v, and possesses additional implicit temperature dependence 
through the first and second derivatives of the pressure and the second derivative 
of the Gibbs' function. Hence, knowledge of these functions is required before 
equation ( 17) can be used for calculating "a" at any state. This point will be 
returned to later. 
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ISENTROPIC TWO- PHASE FLOW THROUGH A HORIZONTAL 
CONVERGING-DIVERGING DUCT 


One of many common equations used to describe the isentropic flow 
of an ideal gas through a horizontal duct is given by 




(18) 


In this equation, A is the cross sectional area of the duct, M is the Mach number 
and V is the gas velocity. With the use of the pseudo-sonic velocity an 
identical equation will be derived for two-phase flow in which a pseudo-Mach 
number will replace its counterpart of ideal gas flow. 


For isentropic two-phase flow of liquid-vapor phases through a 
horizontal converging diverging duct, the continuity equation, general energy 
flow equation and the combined first and second laws are, respectively: 


dA 

A 


dv 

v 


dV 

V 


dh = - 


VdV 

Jg„ 


dh = 


vdp 


(19) 

(20) 

( 21 ) 


Equations (20) and (21) can be combined to give 
VdV 


dp = - 


g c V 


( 22 ) 


Substitution of equation (22) and ds = 0 into equation ( 16) yields 


J<7 


dv = 


g c Tp 


i2 


s av . 

V 


(23) 
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When this is used to eliminate dv in equation ( 19) the result is 



( 24 ) 


This result is then transformed by use of pseudo-sonic relation of equation 
( 17) into 



If a pseudo-Mach number is now defined as 


M = 


V 

a 


(26) 


equation (25) becomes 

dA_ = ( M 2 - 1) dY 
A V 

and this is identical in form to the ideal gas counterpart. 

Equations (27) , (23) , (22) , (21) and (4) yield information of the 
way various quantities increase or decrease depending on whether the flow is for 
M > 1 or M < 1 and further, whether the flow is for a converging section or for 
a diverging section. This is summarized in Table I in which the sign dT 
corresponds to an assumed positive p'^,. 


The analogy with ideal gas flow is carried further in recalling that 
in this flow an improper back pressure on a supersonic nozzle can cause a 
standing normal shock in the diverging section such that the Mach number 
ahead of the shock is greater than one and less than one behind the shock. 
Further the flow is accelerating with falling pressure in the diverging section 
ahead of the shock and diffusing with rising pressure after the shock. Table I 
suggests a similar situation in two-phase flow since both M > 1 and M < 1 
flows, for example, apparently can be accommodated in the diverging section, 
the former characterized by falling pressure and the latter by increasing 
pressure. 
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TABLE 1. EFFECTS OF AREA CHANGE ON VELOCITY, SPECIFIC 
VOLUME, PRESSURE, ENTHALPY AND TEMPERATURE IN 
TWO-PHASE FLOW ABOVE AND BELOW PSEUDO-SONIC VELOCITY 



dA 

dV 

dv 

dp 

dh 

dT 

M < 1 

- 


+ 

- 

- 

- 

M < 1 

4 - 

- 

- 

+ 

4 * 

4 * 

M > 1 

- 

- 

- 

4 * 

4 - 

4 * 

M > 1 

4 - 

4 

4- 

- 

- 

- 


At the throat M = 1; therefore, 

w 

V 2 = a 2 (28) 

l v 

where the t subscript refers to throat conditions. If the mass flow rate per 
unit cross-sectional area is denoted by G the continuity equation then becomes 

V 2 = G 2 v 2 • (29) 

t v L 

Substituting a 2 from equation ( 17) for V 2 in equation (29) gives for the throat 
specific volume after substituting ( 12) for o^: 


v t P t = J *t + 


g p«2 

fc t_ 

G 2 

t 


(30) 


If p" is not zero v^ can be eliminated from equations ( 30) and ( 5) to give a 
throat temperature equation in the form 


7 





No throat subscript is used for the entropy since isentropic conditions are 
assumed. In this form it is again apparent that knowledge of the Gibbs* function 
and the vapor pressure curve are needed to calculate the critical throat tem- 
perature. Further discussions of these functions will be postponed at this point. 


TWO- PHASE FANNO FLOW 


For this adiabatic pipe flow with friction an analogy between two- 
phase equilibrium flow and ideal gas flow will again be developed; in particular, 
it will be shown that choking occurs in the two-phase flow at pseudo-Mach one. 

In this flow the continuity equation is 

V = Gv (32) 


with G now a constant. The energy equation for steady flow is 


h + 


V 2 

2Jg c 


= h 0 


(33) 


where h 0 is a constant. These two equations are then equivalent to 


dh = - 


G* 


J g„ 


vdv . 


(34) 


When equations (34) and (4) are used in the general thermodynamic relation 


ds = ~ dh - — vdp 


(35) 
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to eliminate dh and dp, the result is 


, vp* G 2 v , 

ds = - T7T dT - — — ~ dv 


JT 


J So T 


( 36 ) 


An expression for the differential temperature change is obtained by equating 
the right sides of equations (10) and (36) and rearranging the terms to give 


G 2 v 


+ Tp' 


dT = -I 


J<7 + vp* 

V 


dv . 


(37) 


Elimination of dT in equation (10) by use of equation (37) together with the 
of equations (12), (17) and (26) finally leads to 


use 


ds = 7 — vp ' 2 — 7 [1 - M 2 ] dv . (38) 

J Jcr + vp' 
v 

Since Fanno flow requires a positive ds, equation (38) initiates a sign 
analysis of flows for which M is greater than one and for which M is less than 
one. Table II is a summary of the sign analysis obtained from equations (38) , 
(37) , (34) , (32) and (4) where in the use of (4) , it is again assumed that p'^, 

is positive. Once again the analogy with ideal gas Fanno flow is apparent. 


TABLE II. EFFECTS OF FRICTIONAL LOSSES ON VELOCITY, 
SPECIFIC VOLUME, PRESSURE, ENTHALPY AND TEMPERATURE 
IN TWO-PHASE FLOW ABOVE AND BELOW PSEUDO-SONIC VELOCITY 



dA 

dV 

dv 

dp 

dh 

dT 

M < 1 

0 

+ 

+ 

- 

- 

- 

M > 1 

0 

1 

- 

+ 

+ 

+ 
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Next an integral will be obtained for equation ( 37) that gives the 
temperature, specific volume relation at any point in the flow. This is accom- 
plished indirectly as follows. Equations (32) and (33) are combined to elim- 
inate V, giving 


h + 


G 2 

2Jg c 


= h 0 


and h is then eliminated by the use of ( 1) , resulting in 


(39) 


G 2 -> 

\j/_ + Ts + v* — 

T v 


2Jg„ 


= h. 
— u 


(40) 


Replacing s in equation ( 40) with equation ( 5) then reveals v as a quadratic 
function of T which can be solved explicitly for v, giving at any point in the flow 


v 


g Tp* 
c 

G 2 



2G 2 J (* - T*» - h 0 ) 
g c T 2 p' 2 


(41) 


The right side of this equation is a function of temperature through 
the dependence of p and 4» on the temperature. The choice of sign preceding the 
radical corresponds to the assumption of positive p'^,. 


The entropy and enthalpy dependence on the temperature are then ob- 
tained by eliminating v in equations (5) and (39) by use of equation (41). 


Since M = 1 at the choke point (denoted by the superscript* ) equations 
(32) , (26) , ( 17) and ( 12) can be combined to give 


v* 


gV 


T* 



+ ■ 


(42) 


Then with equation (41) also applied to the choke point, the specific volume may 
be eliminated between equations (41) and (42) leaving one equation to be solved 
for the choke temperature T* and, once again, this solution requires knowledge 
of p,p and . As before, the discussion of these functions will be deferred. 
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TWO- PHASE RAYLEIGH FLOW 


In this horizontal flow of a two-phase mixture only heating effects are 
considered. Fluid friction is neglected and equilibrium is again assumed to 
exist between the phases at all points. Just as in the case of ideal flow, this 
reversible two-phase flow exhibits choking effects at M = 1 and analogous 
characteristics in the two-phase regimes above and below M = 1. 

In this flow a force balance on an element and the continuity equation 
given by equation ( 32) combine to give for the two-phase Rayleigh flow 


P + 



With the use of equation (4) equation (43) becomes 


(43) 


G 2 

p'_ dT + dv = 0 

T «c 


(44) 


and this together with ( 12) substituted into ( 10) gives 


ds = 


1 - 


Jo G 2 ' 
v 

g Tp» 2 
c 


.El 


dv 


(45) 


Equations (32), (17) and (26) then transform equation (45) into 


n | 

ds = (1 - M 2 ) dv 
J 


(46) 


Equations (46) , (44) , (32) , and (4) are summarized with respect to 
sign behavior in Table III. Once again positive behavior of p'^, is assumed and 

the table represents heating of the fluid. For cooling, the signs in the table 
are reversed. 
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TABLE IE. EFFECTS OF HEATING ON ENTROPY, SPECIFIC 
VOLUME, VELOCITY, TEMPERATURE, AND PRESSURE IN TWO- PHASE 
FLOW ABOVE AND BELOW PSEUDO-SONIC VELOCITY 



dA 

ds 

dv 

dV 

dT 

dp 

M < 1 

0 

+ 

4 * 

+ 

- 

- 

M > 1 

0 

+ 

- 


+ 

+ 


For choked flow the asterisk superscript will again be used to re- 
present properties of the fluid at the state M = 1. At this state equation (26) 
reduces to 


f 2 = a* 2 . (47) 

When equations (32) and (17) are substituted respectively for the left and right 
sides of equation (47) the result, after elimination of by equation ( 12) , can 

be arranged to give 


P 


»2 

T* 


+ 



c] 




, = 0 . 
T * 


(48) 


For known functions 


and 4» 

T 


equation (48) can be solved for the 


choke temperature and this then gives the choke pressure and Gibbs’ function, 
p* and 4 e:<. The choke pressure is then used in equation (43) to give the choke 
value of the specific volume. The choke entropy follows from equation ( 5) and 
the enthalpy from equation ( 1). Once the choke pressure and specific volume 
are known an alternate method for calculating the entropy and enthalpy is to 
calculate the quality from the pressure -specific volume values and then use the 
quality to calculate the desired entropy and enthalpy. In any event, equation 
(48) pinpoints the choke point when p and 4' are known functions of temperature. 
In the absence of known equations representing these functions the choke point 
must then be obtained graphically by calculating the Rayleigh states with the 
aid of the usual tabulations of two-phase data. 
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NUMERICAL EXAMPLE OF RAYLEIGH FLOW 


For the Rayleigh flow a numerical example is presented to illustrate 
the use of equation (48) together with the pseudo-sonic velocity and pseudo- 
Mach number concepts when empirical formulations of p and are employed. 

In the example the entrance conditions are saturated water at 400* F. From 
these conditions an iterative procedure is used for the calculation of the 
Rayleigh states in which successively reduced pressures are introduced into 
equation (43) yielding successive values of the specific volume and, hence, a 
series of qualities from which enthalpies and entropies are calculated. These 
results are summarized in Table IV and shown graphically in Figure 1, where 
all of the states except the choke point are calculated in this manner. The 
choke state calculation is initiated with the use of equation (48) and the pseudo- 
Mach numbers are obtained from (12), (17) and (26). For these calculations 
empirical expressions for p and \l> are taken from Bursik and Onuf [1, 2]. 

From Reference 1 the empirical expression for the Gibbs' function for the 
liquid vapor region of H 2 0 is obtained as 

4' = a(r-i) -brlnr (49) 


where 


a = 489. 286 
b = 489.281 


and 


r 



For computational purposes a and b will be taken as equal, i. e. , 


asbs 489. 3 . 


(50) 

(51) 


(52) 


(53) 
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TABLE IV. SUMMARY OF TWO-PHASE RAYLEIGH FLOW EXAMPLE 









































. Btu 

Entropy ~ S lbm . R 

FIGURE 1. DEPENDENCE OF SPECIFIC ENTROPY ON TEMPERATURE 
IN RAYLEIGH TWO- PHASE FLOW 


From equation ( 49) it follows that 




T 


i 

t 3 


[a - b ( 1 + In r)] 


(54) 


or 
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( 55 ) 


*" T = - ITT - • <56) 

It is interesting to note at this point that while a and b are very near- 
ly equal, theoretically they must differ. This can be demonstrated by eval- 
uating both equations ( 54) and ( 55) at the triple point where r = 1. Thus from 
equation ( 54) 


4>' = 

T, 


a - b 


and from equation ( 55) 


*' T = 0 . 


If equation ( 5) is now applied to the saturated liquid state of H 2 0 at the triple 
point where the Keenan and Keyes [3] steam tables have their entropy datum, 
the result is that 


= TV f P* . (59 

X 3 

By the Clapeyron equation p' is greater than zero; hence must likewise 

1 3 *3 

be greater than zero. Equation ( 57) gives a very small positive value for 

which is neglected in equation (58) . 
a 3 

Additional discussion of the Gibbs* function is in Appendix A. 

The summary of Bursik (2) , is included in Appendix B. It shows 
that if some polytropic family of two-phase curves, defined initially by an 
unspecified constant A. as 


A 

pv = K 
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where K is a constant that depends upon X, transforms into the T-s plane as a 
family of parallel curves, then the constant X must be related to a dimension- 
less grouping of thermodynamic variables by the following expression: 




(61) 


The degree to which equation (61) represents two-phase data for the 
liquid-vapor region of H 2 O is discussed further in the Appendix. For present 
purposes equation ( 61) is taken as a sufficiently accurate description of the 
H 2 0 liquid-vapor region when X is given the value 

X = 1.252 (H 2 0) . (62) 

This same polytropic theory yields an additional result that will be useful in 
later numerical examples, namely 


J2. _ h. ~ 1 
P 1 " X 


t 3 


(r - r t ) + 




(63) 


In this result from the polytropic theory the subscript 1 refers to an arbitrarily 
selected reference temperature and r*, pj, v^ ^ and s^^ are all then computed 

at this temperature; hence equation (63) requires p/p' to be a linear function of 
the temperature or the temperature ratio by virtue of (52) . Finally, the poly- 
tropic theory yields an equation for the vapor pressure curve in the form 


T = Tj + 




(64) 


When p" and 4'" are eliminated in equation (48) by the substitution of 
equations (56) and (61) the result is an empirical choke equation 


c + ± _ JG 2 b X 

P T* S c T ® 2 r T* p,2 T* 


(65) 
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At this point equation ( 65) is an implicit function of the choke temperature. 

It can be transformed readily into an explicit function of the choke temperature 
by combining it with equations (63) and (64). However, instead of pursuing 
this path, an alternative cut-and-dry solution was used by employing steam 
tables to solve 


c , . . AG 2 b 1 

" — A + 1 — _ o 

P-p* Jg c ^3 r T* 




( 66 ) 


where p' in equation (65) has been replaced by the Clapeyron relation, 
1 * 

equation (8). 


For the conditions of the example 
G = 536. 5 lbm/sec ft 2 
c = 35780 lbj/ft? . 


Therefore equation ( 66) can be reduced to 


35780 

P T* 


= 2. 252 - 


14. 321 
T* 



(67) 

( 68 ) 


(69) 


As already indicated in Table IV the solution of equation (69) yields a choke 
temperature of 340*F which corresponds to a saturation pressure of 121. 2 psia 
as the choke pressure. 

In the same table the pseudo-Mach numbers were calculated from 



( 70 ) 
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after first calculating v from equation (43). Equation (70) results when equa- 
tions (61) and (56) are combined with equations (12), (17), (26), (8), (61), 
(56) and (31). In conclusion, attention is drawn to the tabulations of the 
pseudo-sonic velocities: a very small value for this quantity at the saturated 
liquid state is indicated; increasing values then accompany falling temperatures 
until a maximum pseudo- sonic velocity is reached; and finally, this quantity 
decreases as the temperature continues to decrease. 

Some of the data of Table IV are plotted in Figures 1, 2, and 3. 


NUMER ICAL EXAMPLE OF NOZZLE FLOW 


In this example of reversible two-phase flow through a converging- 
diverging, horizontal duct, the inlet conditions where a negligible approach 
velocity is assumed are denoted by the subscript 1 and the throat conditions by 
a subscript t. A saturated mixture of H 2 O is assumed to enter at 400* F and a 
specific entropy of 1. 0468 Btu/lbm deg R. This corresponds to an inlet quality 
of 0. 5. 


For this example the velocity of the flow at any section will be 
calculated in two ways: The first is the result of the substitution of steam 
table data into the conservation of energy statement for steady flow, namely 

V 2 = 2Jg (hj - h) . (71) 

c 

The second way will transform the right side of this equation into a 
function of the temperature ratio by incorporating equation (49) and the con- 
stancy of the entropy into it. With the entropy set equal to sj in equation (1) 
and the temperature in the same equation eliminated in terms of r by using 
equation (52) , the enthalpies in equation (71) are replaced by the expressions 
obtained from ( 1) applied to the inlet state and to the variable state without the 
subscript. Finally, equation (49) is substituted into the transformed expres- 
sion for the velocity to give 

V 2 = 2bJg [ (/i + 1 - In r t ) r 4 - (M + 1 - In r) r) (72) 

c 

where m is defined as the constant 
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Btu 
lbm * R 


FIGURE 2. SPECIFIC ENTHALPY AS A FUNCTION OF SPECIFIC 
ENTROPY IN RAYLEIGH TWO-PHASE FLOW 
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FIGURE 3. DEPENDENCE OF THE PSEUDO-SONIC VELOCITY ON 
TEMPERATURE IN RAYLEIGH TWO- PHASE FLOW 
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T 3 S 1 


(73) 


For the assumptions of the example 


491. 69 | 1. 048 
489. 3 


1. 052 


(74) 


_ 859. 69 
ri 491. 69 


1.748 


(75) 


In rj = 0. 5585 . 


With these, equation (72) reduces to 


(76) 


V 2 = 2 (489.3) Jg [2.611 - (2.052 - In r) r ] . (77) 

c 

Equation (71) and (77) are compared in the first two columns of 
Table V, except for the critical point of the throat where equation (72) is 
equated to the square of the pseudo-sonic velocity to locate the critical 
point. This will become apparent when the pseudo-sonic velocity relation 
is specialized to the present example. 


Substitution of equation ( 12) into equation ( 17) gives 


g/p' 2 
a 2 = C 


vp" - J*" 


(78) 


With the entropy held constant at its entrance value s t , equation ( 5) is solved 
for v and the result introduced into equation (78) to give 


Jg (s t + *’) 2 
c 

(Si + 4'') -^7 - 4'" 


(79) 


Then p", 4* ' and 4 f " are eliminated from equation (79) by equation (55), (56) 
and (61) to give 
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TABLE V. SUMMARY OF ISENTROPIC TWO- PHASE 
NOZZLE FLOW EXAMPLE 




376. 18 769 


376 


370 


360 


357. 1 1056 


351. 6 1117 


350 


340 


328 


327. 01 


302. 43 


278 


277. 85 


253. 26 


252 


228. 68 


228 









0. 503 

775. 9 

1.240 

626. 8 

1. 079 

0. 5037 

772.7 

1. 331 

657. 1 

1. 029 

0. 5048 

767.4 

1. 502 

675.7 

1. 001 

0. 5051 

765. 7 

1. 561 

676. 4 ! 

1. 000 



1136 

0. 5060 

762.2 

1.700 

668. 5 

1. 012 

1257 

0. 5066 

756.4 

1. 928 

652. 1 

1. 037 

. 1387 

0.5075 

750. 0 

2.252 

615. 8 

1. 098 



1858 0. 5097 719. 0 4. 549 408. 4 1. 656 



2075 0.5099 702.0 6.829 303.8 


2. 226 



2269 0.5093 685.2 0.233 221.7 3.051 


23 












































































2 




( 80 ) 


Finally, this equation is transformed into an explicit function of the temperature 
ratio by eliminating p/p' from it with equation (63). This then results in 


a 2 = 


bJg (p - In r) 2 

c 

i jj - In r + i. 

A-lr-rj+0 r 


(81) 


In this relation the constant p has already been defined by equation ( 73) and the 
additional symbol 0 represents another constant given by 


0 = 


Piv 


fgl 


X - 1 JT s s 


(82) 


fgl 


A pseudo-Mach number is then formed from equations (72) and (81) , 


giving 


M 2 = 


2[(1+ p-lnrj) r t -(l+ p-lnr) r] ^ _ ln r t i~ 

(p-lnr) 2 |^r-r t +0 X - i r 


(83) 


Substituting the assumed conditions into equation ( 82) evaluates the 
constant as 


1.252 

35613 

1. 8447 

0. 252 J 

778. 26 1 

491.69 1 0.9608 


0. 888 


(84) 


With all of the constants in equation ( 83) now known the pseudo-Mach number 
is set equal to unity for the throat condition and the equation for the throat ratio 
becomes 
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L i. 052 - In r 
3. 968 + — 

r t “ °* 86 r t J 

(85) 

which yields the solution already indicated in Table V as 

t = 357. 1* F . ( 86) 

v 


Since the inlet temperature was 400* F, this represents a temperature drop of 
only 43* needed to establish the throat conditions, considerably smaller than 
that encountered with a diatomic ideal gas expanding from the same tempera- 
ture. 


Some of the results bearing on the accuracy of the empirical re- 
lationships used are shown in Figures 4 and 5. 


DISCUSSION 


The numerical examples used in this paper were not selected for 
any special purpose other than to illustrate the concepts and investigate the 
reliability of the empiricism that was also introduced into these examples. 

Indeed, the empirical relations are necessary before the pseudo-sonic 
velocity and Mach number can be calculated. The empirical calculations 
viewed in the context of exact, steam-table calculations are shown graphically 
in Figures 1,2, and 5. In each of these graphs the choke point calculation 
makes use of the empirical equation for 4^ and p^ while the other points do 

not. Apparently the choke points fit in reasonably well with the other points. 

In Figure 4 a series of points calculated without recourse to the empirical 
formulations is compared to another empirical series. Again the agreement 
between the two appears quite satisfactory. 

In Bursik and Onuf [1, 2] it is shown that the same appropriate form 
of the Gibbs’ function can be used for the liquid- vapor regions of Freon- 12, 
C0 2 , Hg and NH 3 , whereas it fails for sulphur and the solid-vapor regions of 
HgO and C0 2 . 
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O Exact Calculations 
A Empirical Calculations 



Temperature ~ t (°F) 


FIGURE 4. DEPENDENCE OF FLOW VELOCITY ON TEMPERATURE 
IN TWO-PHASE NOZZLE FLOW 


In Bursik [2] the polytropic correlation is shown to give results at 
least as good as the H 2 0 correlation for the liquid-vapor regions of Freon 12, 

Hg, NH 3 , C0 2 , 0 2 and sulphur, and for the solid- vapor regions of H 2 0 and C0 2 . 

The use of the empirical equations for p and 4' permits choke point 
calculations to be made without the use of the quality; hence it is possible to 
obtain solutions outside the physical two-phase range. In the Rayleigh flow, 
for example, the pressure obtained after the choke temperature is calculated 
from equation ( 66) should be entered into equation ( 43) to calculate the specific 
volume at the choke state. This can then be checked against the saturated 
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O Exact Calculation 
A Empirical Calculation 



FIGURE 5. EFFECT OF TEMPERATURE ON THE MASS FLOW 
RATE IN TWO-PHASE NOZZLE FLOW 


specific volumes corresponding to the choke temperature (to determine whether 
a choke solution exists within the two phase region for the assumed inlet con- 
ditions) . 

George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Huntsville, Alabama, August 12, 1966. 
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APPENDIX A. 

AN EMPIRICAL EQUATION FOR THE GIBBS' FUNCTION 


In Bursik and Onuf [1] the empirical approach to an equation re- 
presenting the Gibbs' function begins with the thermodynamic relation 



which is immediately specialized to two-phase regions by recognizing that the 
pressure is a function of the temperature. This permits the partial derivative 
of the pressure to be replaced by p'^ and equation ( 87) then becomes 



j [Tp* - p] 


( 88 ) 


Since the right side of equation ( 88) is a function of temperature only the 
equation may be partially integrated to give 

u = j (Tp' - p) v + f (T) (89) 

where f ( T) is an arbitrary function of temperature resulting from the integration 
and is related to the Gibbs' function as follows. The two-phase constant volume 
specific heat capacity is formed from equation ( 89) by differentiating partially 
with respect to the temperature. This results in 

cr = “ Tp"v + f ' ( T) (90) 

V J 

where primes again denote differentiation with respect to the temperature. Com- 
parison of equation ( 90) with equation ( 12) then provides the connection between 
f and ^ in the form: 

f'(T) = -T¥" . ( 91) 
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In Bursik and Onuf [1] the main emphasis is on seeking correlations 
for the arbitrary temperature function, f (T). Then the Gibbs' function falls 
out as a byproduct through the integration of equation (91). 

By virtue of the Clapeyron relation given by equation ( 8) and the two- 
phase thermodynamic relation 




P v f 
T f g n 


equation (89) becomes 


f (T) - u 



(92) 


(93) 


That the right side of equation ( 93) is indeed a function of the temperature be- 
comes apparent when u and v are replaced by the two-phase relations 


u 



+ x u. 
fg T 


and 


v 



+ x v 


fgn 


These substitutions result in 


(94) 


(95) 


f (T) = u 

T 



(96) 


The right side of equation (96) contains terms that are generally 
tabulated functions of the temperature. Hence the appropriate tables permit 
the graphical study of f (T). An alternate form can be used for graphical 
study by adding and subtracting pv. to the right side to introduce the enthalpy 

*T 

function. This procedure gives 
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(97) 


f(T) = h 

*T 



For the liquid-vapor region of H 2 0, Table VI shows numerical values 
of f computed by use of equation ( 97) in conjunction with the Keenan and Keyes 
steam tables. When f is plotted as a function of the temperature, as shown in 
Figure 6, the points very closely approximate a straight line whose equation is 

f (t) = 0. 9951 t - 31. 848 . (98) 

In this equation the substitution of temperature in Fahrenheit degrees gives f (t) 
in Btu/lbm. The last column in Table VI enables the comparison of the results 
obtained from the empirical equation ( 98) with those obtained from the steam 
table solution of equation ( 97) . 

In terms of the Rankine temperature, equation (98) becomes 

f ( T) = 0.9951 T - a (99) 

and in terms of the ratio defined by equation ( 52) 

f (r) = br - a ( 100) 

where a and b were defined by equations (51) and (52). With these, equation 
( 91) can be transformed to one with r as the independent variable 

d*' k 

+ - = 0 ( 101 ) 

dr r 

which has an integral 

vf* + b la r 3 Cj . (102) 

r 1 

At the triple point r = 1, therefore, the constant is given by 

Cl = (103) 
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TABLE VI. THE FUNCTION f (T) FOR H a O IN THE LIQUID- VAPOR REGION 


t(*F) 

h, 

f = h. - v, 

f V fg f 

f = 

32 

- 0.005 


40 

8. 000 


60 

28.01 


80 

48.00 


100 

67.70 


120 

87.80 


140 

107.7 


160 

127.6 


180 

147.6 


200 

167.5 


220 

187.4 


240 

207.3 


260 

227.2 


280 

247.2 


300 

267. 1 


320 

287.0 


340 

306. 9 



0. 9951 t - 31. 848 

- 0. 005 
7. 956 
27. 86 
47.76 
67. 66 
87. 56 
107. 5 
127.4 
147. 3 
167. 2 
187. 1 
207. 0 
227. 0 
247. 0 
267. 0 
287. 0 
306. 9 
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TABLE VI. (Concluded) 


t(°F) 

h. 

f « h - -£& V 
f V fg 1 

f = 0. 9951 t - 31. 848 

360 

326. 9 

326. 4 

380 

346. 8 

346, 3 

400 

366.7 

366. 2 

420 

386.5 

386. 1 

440 

406. 3 

406.0 

460 

426. 1 

426.0 

480 

445. 8 

445. 8 

500 

465. 5 

465.7 

520 

485.2 

485.6 

540 

504. 8 

505. 1 

560 

524.4 

525. 1 

580 

544. 0 

545. 3 

600 

563. 8 

565.2 

620 

583. 1 

585. 1 

640 

602. i 

605.0 

660 

621. 1 

624.9 

680 

640. 6 

644. 8 

700 

661. 3 

664.7 
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and equation ( 102) becomes 


d* 

+ b In r = 4' , 1 . (104) 

A second integration gives 

+ b(r In r - r) = + c 2 . (105) 

The constant c 2 is evaluated by noting that at the triple point H 2 0 the Gibbs' 
function for the liquid-vapor region is arbitrarily set equal to zero in the 
Keenan and Keyes steam tables by virtue of the fact that both the entropy and 
enthalpy of the saturated liquid are set equal to zero at the triple point tempera- 
ture in these tables. Hence 

c 2 = -b - . ( 106) 


With this result, equation (105) becomes 

= [b + (r - 1) -b r In r . (107) 

r 

It may then be shown that 

b + = a ( 108) 


in the following manner: h^ is first replaced in equation ( 97) by Ts^ and then 

ithe ratio s. /v. is in turn replaced by p'/J by the Clapeyron relation. If the 
fg fg 

resulting equation is written for the triple point temperature where h f is zero 

t 3 

for the Keenan and Keyes tables and then combined with equation ( 59) , the re- 
sult is 


** T = -f(Tg) 


(109) 


which becomes equation (108) in view of equations (99), (100) and (52). Thus 
equation ( 107) is reduced to equation (49) previously used in the numerical 
examples. 


In Bursik and Onuf [ 1 , 2 ] it is shown that this same approach holds approx- 
imately for a number of substances when their liquid-vapor regions are considered. 
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APPENDIX B. 

A POLYTROPIC VAPOR PRESSURE CURVE THEORY 


In Bursik [2] a one- as sumption theory is developed for the pressure- 
temperature equation in any two-phase region by constructing an analog with 
ideal gases. For ideal gases any family of curves defined by assigning a value 
to the polytropic exponent A in the equation 


A 

pv = 


r>.4 


~ X 


( 110 ) 


has the property of transforming into the T-s plane such that each member of 
the family has the same slope at the same temperature as is apparent from the 
slope expression for any member of the A family 


dT 

ds 



( 111 ) 


where c is the well known polytropic specific heat capacity. For a given gas 
A, 

c depends on the heat capacity ratio assumed to be constant and on the identity 
A 

of the family of curves, namely the assigned value of A. For the constant pres- 
sure family of curves c^ becomes c^ and A becomes 0 while the constant volume 

family is characterized by A -*• 00 with c becoming c . In this case as well as 
in the former, it is apparent that the substitution of c^ into equation (111) re- 
quires the constant volume family of curves to have the property that all are 
parallel in the T-s plane at the same temperature. The same is true for the 
constant pressure family and for any other family defined by assigning a fixed 
value to A. 


In the case of two-phase regions the situation is not as simple with 
respect to the transformation into the T-s plane of polytropic families of curves, 
again defined by equation ( 110) . Indeed the T-s plane for H 2 0 as found in the 
Keenan and Keyes [3] steam tables clearly shows that the curves of the constant 
volume family are not parallel to one another at the same temperature. The 
curves of the constant pressure family are parallel in this plane but this is trivial 
since the pressure is a function of the temperature which in turn is one of the 
coordinate axes. 
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If this trivial case is ignored and it is recognized that the polytropic 
family of constant volume curves has lost the parallelism property under 
transformation to the T-s plane, the question is whether there exists one family 
of polytropic curves associated with some one value of X that retains the paral- 
lelism property in the T-s plane. 

To analyze the question the assumption is made that for any two-phase 
region there is one family of polytropic curves that retains the parallelism pro- 
perty under transformation to the T-s plane. When this assumption is imposed 
on two-phase thermodynamics, the polytropic family is further identified as 
the family of curves given by 

vp' = constant ( 112) 

and an expression for the pressure as an explicit function of the temperature is 
developed as is a procedure for calculating X. Various other predictions re- 
sult, all of which can then be tested against tabulated two-phase data. 

The development of the polytropic two-phase theory is initiated by 
rewriting equation ( 10) as 

ds = 7 d(vp') - *"dT (113) 


in which the terms have been regrouped by using the fact that 

d(vp') = vp"dT + p* dv . (114) 


After v is eliminated in equation ( 113) by introducing equation ( 110) the slope 
of any polytropic curve in the T-s plane is given by 


dT 

ds 


X 



1 


_d_ _2L 

dT p 1 /* 


\j»" 


(115) 


In this equation one assigned value of X identifies a family of curves and the 
constant Cj then identifies a particular member of the family. In general, 
assigning a fixed number X does not result in only one value of the slope ex- 
pression at a given temperature even though p and 'I' are temperature functions 
only; since, for these conditions, the number Cj varies from curve to curve of 
the family. The assumption of the poly tropic theory is that from the infinity of 
numbers that can be assigned to X, there is one, namely 
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X = A ( 


( 116 ) 


0 » 

for which the entire family of curves 

p/° = Cl (117) 

transforms into the T-s plane such that the slope of each member of the family 
is the same at the same temperature. From equation ( 115) this can be true 
only if 


_d 

dT IpVAJ 


= 0 


Then the slope of any member of the family of curves will be given by 


dT 

ds 




( 118 ) 


(119) 


Since the Gibbs’ function depends only on the temperature, equation (119) re- 
flects the basic assumption of the theory in that all members of the family 
X = X 0 now obviously have the same slope at the same temperature. Equally 
obvious, the poly tropic family X * X$ must be simultaneously identified as the 
vp* = constant family. At this point, the numerical value of A 0 has not been 
identified. For simplicity the subscript on X will be dropped with the under- 
standing that only the polytropic family of curves X = X 0 is being discussed. 
With this convention equation (118) on the one hand integrates into 


p’ = kp*A 


( 120 ) 


and on the other hand differentiates into the following non-dimensional form 


*4, = X (121) 

PP" 

where both equations hold for all temperatures in the two-phase region con- 
sidered; that is, equation (120) represents an integral of equation (121). If 
Pj and p'i correspond to an arbitrary temperature reference in the two-phase 
region, say Tj, then the constant k in equation (120) can be written as 


38 


(122) 


Pi' 

k = TF 1 

It cannot be evaluated until X is known. Integration of equation ( 120) and evalua- 
tion of the new constant of integration at the arbitrarily selected then gives 

T = T <4 r4 (p^-pi H- 1 ) (123 > 

and this reduces to equation ( 64) when the Clapeyron equation ( 8) is used to 
eliminate p*' in equation ( 122) and k is then eliminated in equation ( 123) . Hence 
the temperature and pressure are explicitly related by either equation ( 123) or 
(64) and all that remains is to evaluate X. This evaluation can be performed 
with the use of equation ( 63) ; therefore this equation will be derived next. 

From equations (120) and (122) the ratio p/p' is given by 



The exponential term is then replaced by equation (64) and pj’ by equation ( 8) , 
giving 


• - -4^- (T - T ‘> + SJf 

fgi 


(125) 


which can be reduced to equation ( 63) by replacing T and T t in terms of r and r t 
from equation ( 52) . The reference state corresponding to the subscript 1 is 
still arbitrary. Replacement of p* in equation (125) by equation (8) yields the 
form amenable to graphical analysis as 


(T-T,) 

S « A 

fg 




(126) 


fgl 


All of the quantities in the ratio on the left side of the equations are tabulated as 


Pv 


a function of T; therefore a plot of the ratio 


Is. 


as a function of the tempera- 


fg 


ture is predicted as a straight line by equation ( 126) . If the prediction holds, 
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the polytropic exponent can be computed after the slope of the straight line 
is determined. 

A second method of determining A is obtained from equation ( 126) 
by substituting for pv^ from the thermodynamic relation 



u. +7 pv, 
fg J fg 


(127) 


When this is done and the result is rearranged, another linear form is obtained 
as 


= A T 

8 fg _X 



Ti 




(128) 


from which A can be obtained graphically. A final graphical determination of 
A can also be made from a log- log plot of equation ( 120 ) with p' replaced by 
equation ( 8 ). 


Figure 7 shows the raw steam table data plotted in the forms 


Hlg. 


and 


Js. 


as functions of temperature. The encircled points are the raw 


fg "~fg 

data points and are shown in relation to the best apparent straight lines drawn 
through them. The data for the plot appear in Table VII. When the method of 
least squares was applied to both sets of data, each slope yielded the same 
value, namely 1. 252. In Bursik [2] a further analysis of these data indicates 
that both sets are curves with inflection points , and an argument is then pre- 
sented showing that the curves theoretically cannot have these inflection points. 
This of course cannot be offered in support of the polytropic theory; but on 
the other hand the data cannot by itself be used to disprove it because of the 
residual data error associated with the experimental inflection points. 
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u f P v f 

TABLE VII. THE RATIOS - St AND FOR H,0 

s, J s» i 

fg fg 

IN THE LIQUID- VAPOR REGION 


t CF) 

PV fg 

1 ? « ,R > 
fg 

7 s - <* R ) 
fg 

50 

26. 85 

482. 9 

100 

33.20 

526. 5 

150 

40. 38 

569. 3 

200 

48. 37 

611. 3 

250 

57. 18 

652. 5 

300 

66.75 

692. 9 

350 

77.00 

732. 7 

400 

87. 86 

771. 9 

450 

99. 19 

810. 6 

500 

110. 8 

849. 0 

550 

122.6 

887. 2 

600 

134. 1 

925. 6 

650 

144.6 

965. 3 

700 

151.2 

1009 
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APPENDIX C. 

THE EFFECT OF INTER-PHASE HEAT TRANSFER 


In single-phase descriptions an infinitesimal disturbance is assumed 
to propagate through the medium in such a way that condensations and rare- 
factions take place so rapidly that there is no opportunity for heat to flow 
from any one element of the fluid to its neighbors. This gives rise to the 
isentropic pressure-density relation in the right side of equation ( 13) . If 
the same argument is used in the case of two-phase flow then an immediate 
time-scale contradiction arises. The constant entropy condition again reflects 
the thermal isolation of any element from neighboring elements but this 
condition simultaneously requires the flow of heat between the two phases with- 
in the element — hence, the terminology, "pseudo-sonic velocity. " 

That the isentropic transformation of a saturated two-phase system 
involves a flow of heat between the phases is illustrated with the use of the 
model in Figure 8. 

In the model it is convenient to think of the closed, saturated, two- 
phase system as being made up of three closed subsystems. In the initial 
state shown on the left, these are indicated by "g" for the vapor phase of 
mass nig and "f" for the liquid phase which is further subdivided into a liquid 

system of mass dm and another mass m^. For the infinitesimally displaced 

system shown on the right it is assumed that the mass dm of liquid has been 
transformed to vapor such that now two "g" subsystems and only one ”f" sub- 
system appear. Only first-order terms are considered for the properties in 
the displaced state. The length dimensions are shown as volumes since they 
are proportional to one another. An energy description for all of the sub- 
systems in the two states is obtained simply by changing the symbol v for u 
wherever v appears in the model. Since W is the total mass of the two-phase 
system the quality is given by 


x 


m 


= — £ 
W 


(129) 


and 


1 - x 


m^ + dm m^ 

M 

w w 


( 130 ) 
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FIGURE 8 . MODEL FOR AN ISENTROPIC TRANSFORMATION 
INVOLVING INTER-PHASE HEAT TRANSFER 

To within first order quantities the reversible first law expressions 
for each of the three closed subsystems become 

dQ g = m g (u g « +ypv g ') dT (13 

dQ fg = (u fg + 7 p V (1J 

dQ f = m, <u.* + 7 pv/) dT 


(13 








where dQ is the heat transferred to the closed vapor system of mass m , 
o g 

dQfg is the heat transferred to the closed system of mass dm that undergoes 
phase change and dQ^ pertains to the closed liquid system of mass ny 


If dm is replaced by dm then the use of equations ( 129) and ( 130) 

g 

enables the set of equations (131) through (135) to be summed into 


dQ + dq. + dQ. = 
g fg * 


+ u. dx + 
fg 


>■*11 

ipdvj 


(1 - x) U.' + X U * 

f g 


dT 


(134) 


The right side of equation ( 135) is recognized as the heat transfer to the com- 
posite two-phase system and is denoted by dQ; that is, 


dQ = dQ g + dQ fg + dQ f 


(135) 


The terms appearing on the right side of equation ( 135) are then 
split into external and internal parts since any subsystem may receive heat 
from the surroundings external to the composite two-phase system as well as 
from an adjacent subsystem within the composite two-phase system. Thus 


dQ = dQ + dQ 
g ge gi 

(136) 

dQ, = dQ. + dQ. . 
fg fge fgi 

(137) 

dQ { - dQ fe + dQ {i 

(138) 


In this manner the formulation approaches that of Prigogine [4] . Here the 
additional subscripts e and i refer respectively to the external and internal 
parts. 


If the composite system is assumed to be insulated, then each of 
the three closed subsystems is likewise insulated from the same external 
surroundings; that is, each of the dQ terms associated with the subscript e is 
zero and for this case equations (136) through (138) reduce to 
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(139) 






(140) 




(141) 


and the right side of equation ( 134) is zero for this constant entropy process of 
the composite two-phase system. Substitution of equations ( 139) , ( 140) and 
( 141) into equation ( 134) results in 


dQ + dQ 

gi 


fgi 


HO 

~ vt fi 


(143) 


While equation ( 142) requires the sum of the internal heat transfers 
to be zero it does not necessarily require the individual terms to be zero. That 
these terms cannot be individually zero is easily demonstrated by assuming that 
they are zero and then showing resulting contradictions. For example, setting 
the right side of equation ( 140) equal to zero requires from equation ( 132) 
that the latent heat of evaporation be likewise equal to zero. Similarly, setting 
equation ( 133) equal to zero results in s^' being zero everywhere on the satu- 


rated liquid curve and equation ( 131) set equal to zero would require s 1 to 


be zero everywhere on the saturated vapor curve. Since none of these results 
is true, the conclusion is that an isentropic transformation of a saturated, 
two-phase system is accompanied by internal heat transfer between the phases 
of the system. The importance of this inter-phase heat transfer will be in- 
vestigated for the simple model of Figure 8 by comparing the expression for 
the pseudo-sonic velocity of equation ( 15) with an expression for the average 
velocity of sound to be derived. 


When Ihe piston in Figure 8 is withdrawn slightly the disturbance 
propagates through the vapor with a velocity denoted by a . The time t for 

6 o 

the disturbance to travel the entire vapor distance L is then 

g 


t = L /a (143) 

g g g 


but 


L /L = x v /v ( 144) 

g g 
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therefore 


(145) 


similarly 


( 1 - x) v. 


where is, of course, the time for the transmuted part of the disturbance to 

travel through the liquid phase. The total time, t, taken for the disturbance 

to travel through both phases is the sum of the two phase times r and r and 

S 

this is given by 


L 

— X 

v 


+ 1 - x) — 

a a 

g f 


If the average velocity of propagation, denoted by a , is taken as the total dis- 
tance traveled divided by the total time of travel, then 


xv (1 - x) v. 


This equation will be compared to the pseudo-sonic velocity equation given 
by a modified form of equation (78) in which is eliminated by equation (56) ; 
p" by (61) and p’ by equation (8) . This results in the following expression for 
the pseudo-sonic velocity: 


V, 
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The comparison of these two equations will be done at a fixed tempera- 
ture of 400" in order to bring out the effect of the quality. At this temperature 
equation (149) reduces to the computational form 


2 = 32.17 (10 6 ) 

22, 428 . 5. 4284 

V \* 


(150) 


Before equation (148) can be reduced to a computational form, the 
pair of numbers a^ and a f in that equation need to be determined. It is assumed 

that these are obtained from equation (150) by respectively replacing v first by 

v and then by v , This, in effect, forces equations (148) and (150) to inter- 
6 1 

sect at x = 0 and x = 1. The justification for this lies in the fact that at each 
of these points only one phase is present, hence there is no interphase heat 
transfer phenomenon. 


With this assumption the computational form of equation ( 148) be- 
comes 


v (10 8 ) 

0. 4283 + 0. 7832 x 


(151) 


and this is compared at a temperature of 400* F to equation ( 150) at eleven 
uniformly spaced values of the quality in Table VUI. Bearing in mind the 
empiricism employed in these computations it appears that the maximum 
deviation between a and a occurs at about a quality of 0. 5. 

As one final note, it emphasized that the a of equation ( 151) corre- 
sponds to an infinitesimal expansion of the original saturated vapor and liquid. 
If the disturbance generated at the piston of Figure 8 is a compression then, 
in the absence of metastable effects, the saturated liquid is displaced infinite- 
simally into the subcooled region and the saturated vapor into the superheat 
region. In this case the saturated sonic velocity corresponding to the single, 
subcooled phase can be about two orders of magnitude higher than the value 
computed at the same point from the two-phase side by use of equation ( 149) . 
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TABLE Vm. COMPARISON OF THE PSEUDO-SONIC VELOCITY 
WITH AN AVERAGE SONIC VELOCITY AT 400* F 


X 

a ft /sec 

a ft/sec 

0 

43.52 

43. 52 

0. 1 

363.6 

400. 9 

0.2 

583.9 

662. 6 

0.3 

758.2 

862.6 

0.4 

905.6 

1020. 0 

0. 5 

1035 

1148.0 

0. 6 

1152 

1253. 0 

0.7 

1259 

1341 

0. 8 

1357 

1417 

0. 9 

1450 

1482 

1.0 

1538 

1538 


49 




REFERENCES 


1. Bursik, J. W. , Onuf, B. O. : Final Report. USARTC Grant DA TC 
44 177 G, June 1963. 

2. Bursik, J. W. : A Thermodynamic Study of Two- Phase Regions. Phd. 
Theses. Rensselaer Polytechnic Institute, June 1961. 

3. Keenan, J. H. and Keyes, F. G. : Thermodynamic Properties of Steam. 
John Wiley and Sons, Inc. , 1936. 

4. Prigogine, I. : Introduction to Thermodynamics of Irreversible Processes. 
Interscience Publishers, 1961. 


50 


NASA- Langley, 1966 M339 


i 


0 * 

I 

I 


" The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
r shall provide for the widest practicable and appropriate dissemination 

I of information concerning its activities and the results thereof ” 

— National Aeronautics and Space Act of 1958 - 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 
important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless 

of importance as a contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: rriformation receiving limited distri- 

bution because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Technical information generated in con- 

nection with a NASA contract or grant and released under NASA auspices. 

TECHNICAL TRANSLATIONS: Information published in a foreign 

language considered to merit NASA distribution in English. 

TECHNICAL REPRINTS: Information derived from NASA activities 

and initially published in the form of journal articles. 

SPECIAL PUBLICATIONS: Information derived from or of value to 

NASA activities but not necessarily reporting the results of individual 
NASA-programmed scientific efforts. Publications include conference 
proceedings, monographs, data compilations, handbooks, sourcebooks, 
and special bibliographies. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 


i 


